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Abstract 

O" 

,—1. ■ A model based on the sudden approximation has been developed to describe 



high energy single nucleon removal reactions. Within this approach, which 

takes as its starting point the formalism of Hansen [1], the nucleon-removal 
Q . 
5 ■ cross section and the full 3-dimensional momentum distributions of the core 

fragments including absorption, diffraction. Coulomb and nuclear-Coulomb 

interference amplitudes, have been calculated. The Coulomb breakup has 

been treated to all orders for the dipole interaction. The model has been 

compared to experimental data for a range of light, neutron-rich psd-shell 

nuclei. Good agreement was found for both the inclusive cross sections and 

momentum distributions. In the case of ^^C, comparison is also made with 

the results of calculations using the transfer-to-the-continuum model. The 

calculated 3-dimensional momentum distributions exhibit longitudinal and 
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transverse momentum components that are strongly coupled by the reaction 
for s-wave states, whilst no such effect is apparent for d-waves. Incomplete 
detection of transverse momenta arising from limited experimental accep- 
tances thus leads to a narrowing of the longitudinal distributions for nuclei 
with significant s-wave valence neutron configurations, as confirmed by the 
data. Asymmetries in the longitudinal momentum distributions attributed to 

diffractive dissociation are also explored. 
PACS number(s): 24.10.-i, 25.60.-t, 25.60.Gc 
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I. INTRODUCTION 

The breakup of light exotic nuclei has proven over the last decade to be a particularly 
useful tool to study the structure of nuclei far from stability (see, for example, [2,3]). In 
the specific case of well developed single-nucleon halo nuclei, such as ^^Be, much insight has 
been gained into the reaction process. Early work by Serber [4], Glauber [5] and Dancoff 
[6] (later refined by Faldt [7]), investigating high-energy deuteron breakup, demonstrated 
that the relevant mechanisms governing the reaction are stripping, diffraction and Coulomb 
dissociation. More recently these concepts have been adapted by Hansen to describe quan- 
titatively the breakup of the one-neutron halo nucleus ^^Be [1,8]. 

The basic premise of the model - dubbed the "sudden approximation" - is that the 
reaction proceeds by the instantaneous removal of a nucleon from the projectile without 
disturbing the remaining nucleons. This approximation is justified for reaction times much 
shorter than the characteristic time for the motion of the nucleons within the projectile. 
A number of assumptions make the calculations particularly simple whilst retaining the 
essential physical concepts: (i) the incident energy is high enough so that the intrinsic 
velocity of the valence nucleon is much smaller than the projectile velocity; (ii) the projectile 
and the fragment follow straight-line trajectories; (iii) final-state interactions are neglected; 
(iv) the tail of the valence nucleon wave function is well developed, so reactions involving 
the valence particle are essentially surface peaked; (v) the target nucleus can be described 
by a "black disk", so that the nucleon scattered or absorbed by the target is not observed; 
(vi) there is only one bound state of the system, the ground state (the completeness of the 
wave functions thus allows the transition probabilities to the continuum to be calculated via 
sum rules); and (vii) only the dominant part (the transverse component) of the momentum 
transfer generated by the Coulomb field of the target is considered. 

These asumptions are well satisfied in the case of reactions involving a well developed 
one-neutron halo nucleus such as ^^Be. In this case, where the rms radius of the halo is 
some three times larger than that of the core, only the asymptotic part of the halo wave 



function was required (see (iv)) [1,9]. This leads to analytical formulae for the transition 
probabilities in impact parameter space and the longitudinal momentum distributions. In 
particular, in evaluating the longitudinal momentum distributions, Hansen proposed that 
to a good approximation the wave function of the valence nucleon could be evaluated at 
the center of the target [9]. This approximation leads to the interesting result that in 
the limit of very small binding energies the breakup cross section factorizes into the free 
neutron-target cross section and the probability that the neutron is at the center of the 
target with a longitudinal momentum close to the incident momentum per nucleon. Thus, 
the model is intimately related to the spectator model of Hussein and McVoy [10]. Whilst 
this approximation is well suited for the evaluation of longitudinal momentum distributions, 
it leads to a strong cutoff in the perpendicular momenta and is, therefore, less appropriate 
for the evaluation of the momentum distributions in a plane perpendicular to the beam 
direction. 

In a recent report, Barranco and Vigezzi [11] have employed a similar approximation and 
obtained the result that the width of the perpendicular momentum distribution essentially 
reflects the target size and is more or less independent of the structure of the projectile 
ground-state wave function. This conclusion does not agree with recent work we have un- 
dertaken whereby the transverse momentum distributions were demonstrated to be sensitive 
to nuclear structure [12,13]. The sudden approximation model requires a proper normal- 
ization of the various transition probabilities, and therefore cannot be extended easily to 
cases where the ground state is dominated hj p or d waves. In such cases, for neutrons the 
asymptotic part of the radial wave function is given by Haenkel functions^, which exhibit 
a strong singularity at the origin and the corresponding breakup probabilities cannot be 
defined properly. Clearly, in such a case, the transition probabilities should be defined in 
terms of realistic wave functions. 



u+ 



hi{iar), where / is the angular momentum and a is the decay constant. 



As discussed in our earlier papers [12,13] and elsewhere [2,14] single-nucleon removal 
reactions appear to be a powerful tool to probe structure far from stability. It is therefore 
highly desirable to test and compare various models of the reaction process in order that 
spectroscopic information may be extracted with confidence. In this spirit we have extended 
the sudden approximation model to deal with single-nucleon removal reactions where the 
ground state has an arbitrary structure for which the single-particle degrees of freedom 
can be disentangled. The model is conceptually simple and almost all observables may be 
calculated with a reasonable numerical effort even when realistic wave functions are used. 

The paper is organized as follows. Sections II to V are devoted to the basic formalisms 
and development of the model. A comparison of the model calculations and experimental 
data for single-neutron removal is described in Section VI. For comparison the predictions 
of the transfer-to-the-continuum model [15,16] for ^^C are also given. The effects of finite 
detection acceptance are discussed in Section VII. Conclusions are drawn and summarized 
in Section VIII. Analytical formulae for various observables are provided in the Appendices. 

II. BASIC FORMALISM 

We assume that the ground state of the projectile {J'^) can be approximated by a su- 
perposition of configurations of the form [/^ ® nljY" , where /^ denotes the core states and 
nlj are the quantum numbers specifying the single-particle wave function of the valence 
nucleon. This is evaluated in a Woods-Saxon potential using the effective separation energy, 
S'^-^-^ = Sn + Ex {Ex being the excitation energy of the core state). Couplings of core states to 
the final state and dynamical excitation of core states in the reaction are neglected. In this 
approximation the reaction can populate a given core state only to the extent that there is 
a nonzero spectroscopic factor C'^S{I^,nlj) in the projectile ground state. When more than 
one configuration contributes to a given core state, the total cross section for single-nucleon 
removal is written, following refs. [17,18], as an incoherent superposition of single particle 
cross sections, 



a_i„(J:) = Y: C'S{I^, nlj) a,,{nlj, S^ff) (1) 

nlj 

The total inclusive single-nucleon removal cross section {(y'^^in) is then the sum over the cross 
sections to all core states. A similar relation holds for the momentum distributions. The 
term agp is the cross section for the removal of a nucleon from a single-particle state with 
total angular momentum j, 

^jm = Rjl{r) Y^ Crailsiyirai{^)XlmM)- (2) 

evaluated with the effective nucleon separation energy S'e// defined above. We consider only 
spin independent transition operators and therefore all formulae are much simpler with the 
wave function, 

M^ = Ri{r)Yi^if-). (3) 

and the normalization, 

" dr\i,om^ = l. (4) 

Certain cross sections may be expressed in a simpler form in a coordinate system traveling 
with the beam. The impact parameter with components in the direction connecting the 
centre of the core with that of the target is denoted by h. Recoil effects are of the order 
O(^), where Ap is the projectile mass number and are neglected here. The situation is 
illustrated in Fig. 1. After the interaction with the target a part of the wave function is 
removed. At this stage the specific form of the wound {{w)) is not important. The removed 
part of the wave function is then, 

^o{x,y,z) if (x,|/,z) g(w) 
6ilj{r) = i 

otherwise. 

The complement (■?/') is defined by the following orthogonal decomposition of the wave 
function, 



ijo = ^ + Sij (5) 

dfilj*6ilj = f dfipSip* = (6) 



If the nucleoli is absorbed, the wave function of the system is locahzed to Sip at the instant 
of colhsion and this is then the state of the remaining core fragment. Consequently, the 
momentum distribution is given by the square of the Fourier transform of Sip. The stripping 
(absorption) probability is given by the volume integral of the wound, 

Paib) = I df]5ij\^ (7) 

As mentioned in ref. [1], "absorption" is defined here in the context of the black disk model. 
The wound complement is normalized as. 



cir-1^|2=l-P, (8) 

The wave function after collision with the target is approximated by, 

^{f)= e"^{^o - 5ij) = e"^i) (9) 

where g is the momentum transfer to the valence nucleon defined by Eq. (A. 4) in Appendix 
A. The sudden transfer of a momentum q attaches a phase e*'^'' to the wave function. The 
physical meaning of the phase is clarified in Section V. This wave function has the normal- 
ization , 

dr%lj*{f)^lj{f) = [ dfi)*{r)i){f) = I - Pa (10) 



Note that ip contains elastic as well as inelastic (breakup) states. The elastic content is given 
by the overlap with the ground state. 



7ei = ] dftProi^ (11) 

The wave function orthogonal to the ground state is, 

ipi{r) =^{r) - 7e;^o(^) (12) 



Now we are ready to construct the wave function for the decaying state (■^/'^(r)), 

M^ = ^o(r) - S^{r) - leie-^'^M^ =4" leic-^'^^o = e-*«^^>i(f) (13) 

As the state which is decaying contains only square integrable functions it is legitimate 
to speak about the norm of this wave function, 

dr\MrW = Jdf]MrW = I - Pa - \lei\^ (14) 

Clearly, all information concerning the absorption of the valence nucleon is contained in the 
wave function of the wound {Sip), whilst ipd will furnish information on the elastic breakup. 
Let us calculate the Fourier transform of the decaying state, 

Mk) = j^, J dre-^'^^Mrl (15) 

The differential cross section in momentum space is given by, 

da 



-^= 27Tbdb\Mkr (16) 

dk -f 

and the total cross section by, 

a = J 27rbdb J dk\Mk)\^ = J 27ibdb{l - Pa{b) - \-fei{b)\^) (17) 

Now it is clear that the quantity, 

Pelib) = 1 - Pa{b) - \jel{b)\^ (18) 

should be interpreted as the probability in the impact parameter representation of the elas- 
tic breakup process. Relation (18) shows that the sudden approximation model accounts 
for stripping (Pa), elastic breakup or dissociation (Pei) and nuclear and Coulomb elastic 
scattering {'~fei)- Other inelastic processes such as core absorption, simultaneous absorption 
of the valence neutron and the core [19] or "damped breakup" [20] are not included in this 
model. We also note that the closure relation (18) illustrates the need to use realistic wave 
functions in order to obtain properly normalized breakup probabilities. 



III. BREAKUP PROBABILITIES 

In this section we discuss in detail the calculation of the breakup probabilities using two 
specific forms for the wound {w). Using the same notations as in ref. [1] the planar cutoff is 
defined as, 

(w) : (x, y, z) e R^, X > bi = b - Rt 

The second, the so called cylindrical wound is defined for arbitrary z as, 

{w) : (x, y, z) e R^ (x - b)2 + y' < Rt 

where Rt denotes the target radius. For clarity it is useful to define the (normalized) valence 
density by averaging over projections of the angular momentum m, 

PUr) = IT.\H" = r^iir) (19) 

Writing f = {x,y,z) = {s,z), we define also one and two-dimensional projections of this 
density as, 

dydzpyai = 271 rdrp^aiijr) (20) 

p{s) = J dzp^aiir) = 2^ -7^f^^Pvai{r) (21) 

The singularity in (21) is weak and can be integrated by parts. Without loss of generality 
one can assume refiection symmetry, Pvaii^f^) = Pvaiij'), then the mapping x -^ p{x) is a 
refiection symmetric homomorphism, p{—x) = p{x) and 2/q°^ dxp{x) = 1. Since the valence 
density p^i is normalized to one, the density p also satisfies the closure relation / dsp{s) = 1. 

A. Elastic probability 

According to Eq. (11) the elastic probability defined by, 

7d = / dfroi^e^'^i^o - H) =lc- Ic+N (22) 



decomposes into a Coulomb and a nuclear + Coulomb amplitude. It is easy to show that, 



'-fc = '2' dx cos qxp{x) (23) 

Jo 

^c+N = 27r / dxe'^^^Pix) (24) 



For a cylindrical wound the Coulomb and nuclear amplitude is somewhat more complicated 

7C+7V = / ' sds r #e*'?(^+^=°^'^)p(y'62 + ^2 + 26s cos0) (25) 

Note that 7c+Af is a complex function and depends on the specific form of the wound, 
whilst the Coulomb amphtude is real and wound independent. Moreover limq_>o7c = 1 
which shows that in the absence of a Coulomb field or at very large impact parameters 
the Coulomb amplitude is unity. However, the long range Coulomb interaction make this 
convergence very slow. In the case of a well developed neutron halo, the wave function varies 
little over the wound region and one can replace in (25) p by some average value po which 
leads to, 

bc^N? = pt^Jl{qRt) (26) 

in analogy with the Fraunhofer diffraction pattern produced by an absorbing disk. This 
approximation illustrates the diffraction content of 7 [11]. In practice this amplitude is 
evaluated using the exact equations (24) and (25). 

B. Absorption/stripping probability 

The stripping probability is evaluated using Eq. (7), 

Pa{h) = / dxdydzp^ai{r) (27) 

J{w) 

which leads after some simple manipulations to, 

/■oo 

Paib) = / dxp{x) (28) 

Jbi 

in the case of a planar cutoff and , 
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/■Rt /•27r I 

Pa{b)= sds rf0p(J62 + s2 + 2focos0) (29) 

Jo Jo * 

for a cylindrical wound. Independent of the form of the wound we have Pa{0) = 1/2 and 

in the limit g ^ 0, •yc+Nib) = Pa{b) and Pei{b) = Pa(b) — P^ip). The total cross sections 

for stripping and diffraction are obtained by integration of the above probabilities over the 

impact parameter with the volume element 2'Khdh. The above relations imply that for a light 

target, where the Coulomb component is negligibly small, agtr ~ (^diff- In general, however 

astr > o'diff in agreement with the original formulation of Glauber [5] . 

IV. MOMENTUM DISTRIBUTIONS 

As final state interactions are neglected, the momentum distributions in the coordinate 
system traveling with the beam are given by the square of the Fourier transform (Eq. (15)) 
of the wave function Eq. (13). The three components in Eq. (13) lead to an amplitude of 
the form, Aim{k) = Aq + As + Ac+n, where Aq is the unperturbed (intrinsic) amplitude. 
Using the same notations as in ref. [1] one has for angular momentum (Im), 

"^^M^k) = -n^:^i'YUk) j^ r^drji{kr)Ri{r) (30) 

and 

Ac+NMnik) = -lei{b) .^ j Yimik,) / r^drji{kgr)Ri{r) (31) 

[Zti)'^'^ Jo 

— * 

where kg = {kx + q,ky,kz)- The calculation of the amplitude As is more involved. The 
difficulty arises from the angular part of the wave function. The simplest way to proceed is 
to express the spherical harmonics in Cartesian coordinates (see for example ref. [21]). For 
the planar cutoff case we need to evaluate integrals of the form, 

_, /"OO /"OO /"OO 

Ipstik) = / dxxPe-'''^'' / dyy'e''^yy / dzz'e-'''^'Ri{r)/r^ (32) 

Jbi J—oo J—oo 

where p, s, t are positive integers. In the simplest case, p = s = t = 0, we have, 

Jooo(A;) = 27r / rfxe"^'^-^ / uduJo{k±u)Ri{Vx^ + u^r^ (33) 

Jbi Jo 
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where Jq is the cyhndrical Bessel function and A;_l = Jky + kl- For any other combination of 
p, s, t, the corresponding integral is obtained by the appropriate parametric differentiation 
with respect to kx,ky,kz. For an / = wave function the result is, 

^As{k) = 1= / dxcoskxX / uduJo{k±u)Ro{r) (34) 

V2(27r)2 Jbi Jo 

'!^As{k) = —^ / dxsinkxX / uduJo{k±u)Ro{r) (35) 

V2(27r)2 J6i JO 



with r = v^?+a?. The corresponding expressions for p (/ = 1) and d {I = 2) states are 
more involved and are detailed in Appendix B. 

Once the amplitudes have been computed, the differential cross section is obtained by 
averaging over the magnetic projections and the impact parameter, 

% = iT.fdb\AUh,ky,k.)\' (36) 

CtrC L m 

It can be seen that the reaction mechanism modifies substantially the momentum content 
selected by the reaction. For example, for / = the unperturbed amplitude Eq. (30) is 
spherical and real and the amplitude selected by stripping is asymmetric and complex. 
These effects will be discussed in detail in the next section. It should be noted that all 
amplitudes tend to zero for 6 -^ oo as is evident from Eqs. (34) and (35). Furthermore 
the amplitudes (30) and (31) become identical for large h as in the limit 6 ^ oo, 7ej = 1, 
kg = k and the two amplitudes essentially cancel each other. Asymptotically, the amplitude 
for diffraction is, therefore, essentially the Fourier transform of the wound. Both stripping 
and diffraction amplitudes behave for large kx as 1/kx. The reason for this is of course that 
the leading term of the Fourier transform of the step function (characteristic of the black 
disk approximation for the S-matrix ) behaves like l/kx- This property is independent of 
the asymptotic behavior of the wave function and merely characterizes the sharp edge of the 
target. 

In more realistic reaction models, such as the Glauber model with the S-matrix generated 
by an optical potential [20,12,13], the absorption evolves smoothly from to 1 and the 
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amplitudes fall off faster. One may also note that the differential cross section (36) is 
calculated in the intrinsic reference system for a particular position of the target. In this 
system the diffraction amplitude is not symmetric with respect to k^ since the Coulomb field 
pushes the core in one direction. The amplitudes, however, are symmetric to rotations in the 
(x, y) plane. This is obviously not an observable symmetry. In order to obtain observable 
{p\\,p±) distributions one should average over all directions in the {x,y) plane as illustrated 
in Fig. 2. In practical terms, the physical momentum distributions are obtained from the 
following unitary transformation. 



-— = —- df3^{px cos (3 - py sin (3, p^; sin (3 + Py cos /3, p^) , (37) 

dp In J dk 



where the angle f3 is defined in Fig. 2. Equation (37) gives the 3-dimensional momentum 
distribution in the laboratory system. It contains all the physics of the process and together 
with (36) provide a means to study the effects of finite detection acceptances and other ex- 
perimental effects on the momentum distributions. Unfortunately the corresponding calcu- 
lations are very time consuming. In practice we have combined competitive Gauss-Legendre 
numerical integration (up to order 98 in the approximating polynomial) with Monte Carlo 
simulations in the last step (Eq. (37)) which took into account all experimental broadening 
effects: angular dispersion in the secondary beam, angular and energy straggling in the 
target and detector resolution [13]. The Lorentz broadening is also included, 

EB = -f{EA-0-pA), (38) 

PB=PA + 7/3(^^-i5A), (39) 

7 + 1 

where the reference system A is traveling with the beam velocity (3 and B is the laboratory 
system. 

If the acceptance in the plane perpendicular to the beam direction is infinite, then it 
is possible to obtain much simpler formulae for the longitudinal momentum distributions. 
These are detailed in Appendix C. 
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V. COULOMB DISSOCIATION IN THE SUDDEN APPROXIMATION 

As already remarked in ref. [1], the Coulomb excitation calculation based on the (non- 
perturbative) sudden approximation is consistent with perturbation theory [22] in the sense 
that it contains its leading term as a limit valid for low Z^Zt (core and target charge) and 
high velocities. This approximation is no longer valid for large Z^Zt and is thus limited to 
low Z targets. Moreover, in this model the dependence on the projectile energy is given by 
the Coulomb contribution alone via the momentum transfer q. In more sophisticated models 
of the Glauber type (see, for example, [13] and references therein), there is an additional 
energy dependence through the S-matrix elements of the nucleon-target and core-target 
interaction. This dependence, however, is expected to be weak in a high energy regime. 

In this section we wish to clarify the Coulomb dissociation effects present in the model. 
In this respect we make contact with the more elaborate calculations of ref. [23]. This 
model solves the time dependent three-body problem assuming that the core moves along 
a classical (straight line) trajectory R{t) = & + vtz^ whilst the valence nucleon is subject to 
the interaction, 

V^2(r, t) = VUr + m) + Vd^p{f, t), (40) 

where the first term is the neutron-target interaction and the second is a dipole approxi- 
mation for the "shake-off interaction" . If only the Coulomb shake-off is taken into account 
then for a valence neutron. 

In the general case one uses the effective dipole charge defined in Appendix A. Note that 
Eq. (41) embodies a pure recoil effect (~ 1/^p) which retains both the transverse and 
longitudinal momentum transfers to the neutron arising from the projectile-target Coulomb 
interaction. 

The transfer-to-the-continuum model breakup amplitude takes the form. 
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gUKh) = - dt< 0^m(r,t)|0L > (42) 

til J —oo 

where (p^^^^ are initial (final) state wave functions and Im are, as before, single-particle 
quantum numbers carried by the ground state. In the eikonal approximation, we have 

9Uk,b) = ^|rfrYrfte-^'^^+^-*e*r'^*'^^(-^'*V2(f,t)0,„(f) =< k\I{u)\lm > (43) 

where I{uj) is the time integral and huj = e^ — eo is the excitation energy. Pure Coulomb 
effects in the dipole approximation are obtained by swiching off the neutron-target nuclear 
interaction Vnt, 

Icioo) = ^r dte'-'e^ r '^*'^--(*V,,p(t). (44) 

th J-oo 

In the sudden approximation (sa) the excitation energy is negligibly small {uj ~ 0). The 
time integration is thus trivial and we find, 

jsa ^ ^-ixc _ I (45) 

with, 

1 /-"^ ,,, ,, I f^ , ZcZte^ sb + zvt 2ZcZte'^ sb 

where the neutron position vector relative to the core has been decomposed into transverse 
and longitudinal components (s, -z). Comparison with Eq. (A4) in Appendix A shows that, 

e'""^ = e'"^' (47) 

and the complete equivalence of our amplitude Aq+n (Eq- 31) with the corresponding sudden 
approximation amplitude in the eikonal TC model. Therefore, the amplitude Ac+n contains 
breakup to all orders in the dipole shake-off Coulomb interaction at zero excitation energy. 
It is appropriate for evaluation of Coulomb effects for light targets. For heavy targets, the 
radial integral in Eq. (31) is difficult to evaluate numerically. We further stress that the third 
term in the wave function (Eq. (13)) 7e/e*'^^-?/'o is nothing other than the high-energy eikonal 
approximation to the scattering wave function where in this particular case the eikonal 
phase due to the dipole (shake-off) Coulomb interaction is calculated along the unperturbed 
(straight line) classical trajectory. 
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VI. DISCUSSION 

As detailed in refs. [12,13] we have previously undertaken an experimental study of high- 
energy one-neutron removal reactions on a range of psd-shell nuclei. In the present section 
the inclusive cross sections, longitudinal and transverse momentum distributions obtained 
for reactions on a carbon target are compared to the results of calculations using the model 
developed in the preceeding sections. 

As in our earlier work [12,13] the spectroscopic amplitudes entering in Eq. (1) have 
been calculated with the aid of the shell- model code OXBASH [24]. Where known, the 
experimentally established spin-parity {J'^) assignments and core excitation energies have 
been used. In all other cases shell model predictions were employed. The carbon target 
radius was fixed to Rt = 1.15(12)^/^. The core radii have been evaluated with a liquid drop 
formula [25], 

with Xi = 1.17, X2 = 1.225 and 0:3 = —0.115. The single-particle wave functions were 
obtained by solving the Schrodinger equation for a Woods-Saxon (WS) potential including 
central and spin-orbit terms. The depth of the central potential was adjusted to reproduce 
the known effective neutron binding energy (5'^'^'^ = Sn+E^). The potential radius was taken 
to be equal to the core radius and the diffusivity was fixed to a^s = 0.6 fm^. The minimum 
impact parameter was defined by fe^m = Re + Rt and the maximum impact parameter 
was fixed to be bmax = 50 fm which ensured that there were no spurious effects from the 
cancellation of the amplitudes of Eqs. (30-31). Using the upper adiabatic limit suggested by 
Hansen [1] produced only minor changes in the cross sections. We further assumed that the 
core ground and excited states have the same density distributions and, as such, the same 
Woods-Saxon geometry was employed for all core states of the projectile. 



^A fine tuning in the range 0.5 — 0.65 fni could, in principle, be made for each nucleus to improve 
the cross section with respect to the measured value. 

16 



The cross sections calculated within the planar cutoff approximation are displayed in 
Fig. 3. The stripping and diffraction (including Coulomb brealcup) components are also pre- 
sented. A good overall agreement with the experimental cross sections is observed. Keeping 
all the parameters fixed as defined above, the results obtained with a cylindrical wound are 
shown in Fig. 4. This calculation systematically underestimates the experimental value by 
a factor of two. This is a pure geometrical effect resulting essentially from the fact that the 
cylindrical wound underestimates the real extension of the interaction region. As such a 
diffuse edge cylindrical wound would probably be more appropriate. 

The particular form of the wound, however, has essentially no influence on the shape of 
momentum distributions. Calculations in the planar cutoff approximation are displayed in 
Fig. 5. Monte Carlo filtering of the calculated distributions including Lorentz broadening, 
energy straggling in the target, the emittance of the beam and detector resolutions have 
been performed and the resulting distributions, normalized to the peaks of distributions are 
displayed by the solid lines in Fig. 5. In addition to the excellent overall agreement, we note 
that the distributions are somewhat better reproduced by the sudden approximation model 
as compared to the Glauber calculations employed in our earlier work [12,13], as illustrated 
in Figs 6 and 7. This is primarly attributable to the strong damping of the high momentum 
components in the wave function by the sharp cutoff approximation. 

As mentioned earlier, the momentum distributions are first calculated in the projectile 
rest frame which is determined by a particular relative core-target configuration. The observ- 
able momentum distributions are then obtained by the unitary transformation of Eq. (37). 
The effect of this transformation is displayed in Fig. 8 for the the s-wave valence neutron 
in the ground state of ""^^C (S„=1.2 MeV). In the rest frame, the kx and ky distributions 
have markedly different shapes as the Coulomb shake-off imparts momentum in only one 
direction (x) . After averaging over all directions (Eq. 37) the Px and Py distributions become 
identical in agreement with what is observed experimentally. The asymmetry induced by 
the Coulomb interaction in the rest frame, translates into a small broadening effect in the 
laboratory frame. 
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Further insight into the role played by the reaction mechanism on the momentum distri- 
bution is explored in Fig. 9. Here momentum distributions in the transparent limit of the 
Serber model [4] are compared with the planar cutoff calculations in the rest frame for s and 
(i-wave valence neutrons. As already noted by Hansen [9], the longitudinal components (kz) 
become narrower, irrespective of the angular momentum carried by the wave function. For 
the transverse component (kx), the effect is more complex. A broadening effect is observed 
for the s-state, whilst the shape of the distribution is completely changed for a rf-state. The 
overall effect after averaging (Eq. 37) is a broadening of the transverse distributions and a 
narrowing of the longitudinal one. These conclusions are consistent with the measurements 
of ref. [13] whereby the transverse distributions were found to be sytematically broader than 
the longitudinal ones. We note that this is in contradiction to the analysis of Sagawa and 
Takigawa, whereby (for s-wave states) the transverse distribution is narrowed by the ab- 
sorptive cutoff induced by the reaction process [26]. The calculated transverse momentum 
distributions, after inclusion of the experimental effects, are compared to the data for se- 
lected nuclei in Fig. 10'^. The finite angular acceptances of the spectrometer and detector 
resolution introduce a smooth cutoff in the high momentum tails of the distributions. As 
in the case of the longitudinal momentum distributions, the shape and width of transverse 
distributions show a direct dependence on the projectile structure, and are not simply a 
reflection of the target size as suggested by Barranco and Vigezzi [11]. 

As a further example we display in Fig. 11 the results of calculations for ^^C for the 
three possible ground-state spin-parity assignments. The results are very similar to those 
obtained in our earlier work using an extended Glauber-type model (Fig. 19 of ref. [13]), 
which suggested a spin-parity assignment of 3/2+ arising from a dominantly ^^C(2j'" ® i'dr^/2) 
configuration. As noted in refs [12,13], this is in line with the direct observations of coincident 



^Owing to the very time consuming nature of these calculations, only px distributions for four nu- 
clei with ground state structures representative of those measured in ref. [13] have been computed. 



1.76 MeV 7-rays by Maddalena et al. [27] and more recently by Datta Pramanik et al. [28]. 



For completness, the perpendicular momentum distributions (p_L = Jp^. + p^) have also 
been reconstructed [32] from the data of refs. [12,13]. The calculations for selected nuclei 
are compared, after Monte Carlo filtering of the experimental effects to the data in Fig. 12 
where very good agreement is again found. 

As intimated earlier, one of the goals of the present work was to explore a reaction 
model other than the Glauber-type approach which has been the principle means to utilising 
high-energy nucleon removal as a spectroscopic tool. In this spirit we have also performed 
calculations for the reaction of ^^C on a carbon target using the Transfer-to-the-Continuum 
Model (TCM) developed by Bonaccorso and Brink [15,16,29], which has been employed to 
describe with some success single-neutron removal from beams of ^^'^^Si and ^^S [30]. As 
described in ref. [29], the TCM formalism includes a more complete dynamical treatment 
of the motion of the removed nucleon than Glauber/eikonal models. In particular, spin- 
coupling between the initial and final states may be included and can result in asymmetric 
momentum distributions. 

The TCM calculations presented here (Fig. 13) were performed as outlined in ref. [29] 
(Eqs. (2.1-2.4)). The neutron-target optical potential was taken from ref. [31] and the strong 
absorption radius, used in the parametrization of the core survival probability, was fixed to 
be 6.8 fm. The ^^C core longitudinal momentum distribution was obtained from that of the 
neutron by momentum conservation. As may be seen in Fig. 13, the results assuming the 
^^C ground-state structure outlined above (J^^=3/2"'") are in excellent agreement with both 
the data and the calculations using the sudden approximation. 

VII. DETECTION ACCEPTANCE EFFECTS 

Following the first measurements of core-fragment longitudinal-momentum distributions 
[33], the infiuence on the observed lineshapes of limited angular or transverse momentum 
acceptances has been discussed [34-38]. These discussions have, however, been based on the 
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assumption first introduced by Riisager of a 3-dimensional Lorenzian momentum distribution 
[34]. More recently it has been conjectured that the reaction mechanism itself causes the 
core-fragment momentum components to decouple [9]. As such, incomplete detection in the 
plane {x, y) perpendicular to the beam direction would have no infiuence on the measured 
longitudinal momentum distribution. 

The data set presented in refs. [12,13] provides a good opportunity to investigate such 
an effect at a quantitative level, as the full 3-dimensional momentum distribution of the 
core-fragment have been acquired. In the off-line analysis, the angular acceptance of the 
spectrometer was reduced from the full acceptance of ±2° to ±1.5, 1, 0.5, 0.25°. These limits 
correspond to transverse momenta p^ and Py of around ±200, 150, 100, 50 and 25 MeV/c. 
Of the 23 nuclei studied in refs. [12,13] the restricted acceptances lead to a narrowing of 
the longitudinal momentum distributions in 5 cases — ^^B, i5,i6q g^j^^^ 24,25p^ jj^ ^j^g ^^^^ 
of ^^C the narrowing reaches some 25 % for the most limited acceptance (Fig. 14 and 
15). Interestingly these nuclei have in common a large s-wave component in the ground- 
state wave function. For other nuclei, with ground states dominated by rf-wave valence 
neutron configurations, no significant reduction in the widths of the longitudinal momentum 
distribution was observed. 

In the sudden approximation model the full 3-dimensional core-fragment momentum dis- 
tribution can be evaluated and projected onto the desired direction after integrating over the 
corresponding acceptances. For the four nuclei selected earlier as examples, the acceptances 
applied to the data have also been applied to the calculated distributions. The narrowing 
with the reduced transverse acceptances of the longitudinal momentum distributions is well 
reproduced for ^^B and ^^'^^c (Fig. 14). In the case of ^^N, experimentally the width di- 
minishes by only some 5% as the acceptances are reduced. This trend is well reproduced 
by the calculated distributions^. We note that, as intimated above, the ^^N ground state is 



^Note that the axis displaying the FWHM in Fig. 14 has been expanded and the widths are in 
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dominated by a rf-wave valence neutron configuration. 

These calculations indicate that the momentum components are, contrary to the sugges- 
tion of Hansen [9] strongly correlated and reduced acceptances in, for example, the transverse 
direction will affect the longitudinal momentum distribution. This effect is most evident for 
s-wave valence neutron configurations (e.g. ^^B and ^^'^^C). The longitudinal momentum 
distributions for each acceptance are compared in Fig. 15 for the case of ^^C with the cal- 
culated distributions and very good agreement is found. 

For the distributions derived using the full acceptances of the spectrometer, the experi- 
mental distributions are somewhat asymmetric and exhibit low momentum tails. This effect 
is not reproduced by the present model and, as noted in our earlier paper [13], almost 
certainly arises in the case of weakly bound systems as a result of a strong coupling to con- 
tinuum in diffractive dissociation. Such a process has been succesfully modelled by Tostevin 
et al. [39] within the coupled discretized continuum chanels (CCDC) formalism, where the 
redistribution of relative energy to the internal excitation of the projectile is treated exactly. 
As is evident from Fig. 15, as the angular acceptances are progressively reduced the asym- 
metry becomes less pronounced as the contribution from diffraction decreases, an effect also 
observed by Tostevin et al. [39]. 

Finally, the case of ^^N is displayed in Fig. 16. As noted above, within the statistical 
precision of the present measurements the width of the longitudinal momentum distribution 
remains unchanged with decreasing acceptances. This supports the model prediction that 
for (i-states the different momentum components are effectively decoupled. Interestingly the 
asymmetry in the momentum distribution appears to persists even for very limited transverse 
acceptances a feature which cannot be easily explained by diffractive processes. 



fact reproduced to within 5% or better 
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VIII. CONCLUSIONS 

The sudden approximation approach for the description of high-energy, single-neutron 
breakup of halo nuclei, originally introduced by Hansen [1], has been extended to include 
realistic wave functions and to incorporate shell model spectroscopic amplitudes. The theory 
is based on the strong absorption description of the core-target and neutron-target inter- 
actions. Applied to single-neutron removal reactions, the model allows for the calculation 
of the full 3-dimensional momentum distribution of the core-fragment including stripping, 
diffraction and Coulomb dissociation mechanisms. The Coulomb/nuclear interference is also 
taken into account. As in the case of our Glauber-type calculations [12,13] the observation 
that the transverse momentum distributions are systematically broader than the longitudinal 
distributions is reproduced by the present calculations. 

Calculations were performed for comparison with measurements of inclusive cross sections 
and longitudinal and transverse momentum distributions for a series of some 23 neutron- 
rich p — sd shell nuclei [12,13]. Suprisingly, for such a relatively simple model, very good 
agreement with the measured cross sections and momentum distributions was found. Indeed, 
the momentum distributions were somewhat better reproduced by the present model than 
the more sophisticated Glauber-type calculations [12,13]. 

Importantly, the effect of limited detection angular acceptances on the longitudinal mo- 
mentum distributions could be investigated. A significant reduction of the widths was 
observed experimentally for nuclei with ground states dominated by s-wave valence neutron 
configurations. Little or no reduction was observed, however, for nuclei with dominant d- 
wave valence neutron components. This effect was well reproduced by the present model 
calculations and is believed to arise as a consequence of the correlations between the momen- 
tum components in the 3-dimensional momentum distribution of the core-fragment following 
the reaction. 
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APPENDIX A 



In this appendix we derive a convenient expression for the classical momentum transfer 
using straight line trajectories and retaining only the dipole component of the Coulomb field. 
The projectile consists of core {rric, Zc) and a cluster {rrix, Z^) moving in the 2;-direction with 
velocity v. The impact parameter b with respect to the target is measured in the x-direction. 
The dipole effective charge is defined by [40] as, 



7(1) 
'eff 



rric + m^ 



(Al) 



The time dependent electric field due to the target charge Zt is given by, 

iZte 



U^ 



E{t) 



(52 + ^2^2^2)3/2 







yvtj 



here 7 is the Lorentz contraction factor. The classical momentum transfer is then. 



Ap 



dteZiffE{t) 



(A2) 



One sees immediately that Ey = ^ Apy = and Ap^ = from parity consideration. It 
follows that the momentum transfer has only one component in the x-direction given by, 

dt 



Ap, = ^ZtZi^e'b 



2Z,zi]\e' 



(62 + ^2^2^2)3/2 
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(A3) 



The momentum transfer imparted by the cluster is , 



(A4) 



2ZtZ^^}y b 



hep 62 

Momentum conservation gives then the core momentum. When the cluster is a neutron 
Eq. (A4) is identical to Eq. (10) of ref. [1]. 

APPENDIX B 

Explicit expressions for As amplitudes (Section IV) are detailed in this appendix, for 



I =0, 1 and 2 wave functions. In all expressions the notation r = VM^ + a?, is used. 
For / = one has, 
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For / = 1 and m = ±1, 
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For I = 1 m = 0, 
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Finally, for I = 2 m = one has, 



with 



A. 



<520 



Stt 



(3/i + h) 



1 k^r^ f i^fty 3 

r"o~ / QjX cos rC^x 







e 



1^^ roo rmtty 

^Ii = — yy / dxcoskxX / u^duJi{k±u)R2{r 
u^du{Jo{k±u) — J2ik±u))R2{r 

S/i = -r4- dxsinkxX / u^duJi(k±u)R2{r 
k'l^ Jbi Jo 

1 k'^ roo I- infty 

-\ — -^ / dxsinkxX / u du(Jo(k±u) — J2(k±u))R2(r 

2 kj_ Jbi Jo 

/•OO /"OO 

^h = / dxcoskxX / u^duJo{k±u)R2{r 

Jbi Jo 

/•OO /"OO 

Q'/2 = — / (ixsin/cj-x / u^duJo{k±u)R2{r 

Jbi Jo 



(B21) 



)/r^ 




)/r2 


(B22) 


)/r^ 




)/r^ 


(B23) 


)/r^ 


(B24) 


)/r^ 


(B25) 



APPENDIX C 



In this appendix we give explicit formulae for the longitudinal momentum distributions in 
the case of a planar cutoff. Complete transverse acceptances are assumed. The momentum 
distribution is given by. 
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where $ is a generic notation for Sip in the case of stripping or ipd in the case of diffraction. 
The differential cross section is then, 
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1 = 1 m = ±1 
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Note that the apparent singularity in the above integrals is removed by the behaviour of 
the wave function near origin Ri{r) ~ r'+^. For diffraction the same expressions hold except 
that the integral over x is replaced by, 



dx\f{x,b)\' 



(C9) 



where 



f{x,b) = l-9{x-b,)-jei{b)e- 



iqx 
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6i,7e/ and q are defined in the main text. ^ is a step function. -4 is already included 
in above expressions. The total differential cross section is obtained simply by summing 
up contributions from various m values. For a given / state, the longitudinal momentum 
distribution is symmetric. 
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FIG. 1. Schematic view of the "wound" induced by the projectile-target interaction for (a) a 
planar cutoff and (b) a cylindrical wound. 




FIG. 2. Schematic representation of the unitary transformation of Eq. (37) of the momentum 
vector k in the rest frame to the laboratory momentum (p) by rotation through an angle /3. 
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FIG. 3. Single-neutron removal cross sections calculated in the planar cutoff approximation 
(open circles, dashed line) compared to the experimental values [12,13] (full circles, solid line) for 
reactions on a carbon target. The contributions arising from stripping (open triangles, dash-dotted 
line) and diffraction plus Coulomb dissociation (open diamonds, dotted line) are also detailed. 
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FIG. 4. The same as in Fig. 3 but for a cylindrical wound. 
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FIG. 5. Core- fragment longitudinal momentum distributions for reactions on a carbon target. 
The experimental data are taken from ref. [12,13]. The solid lines correspond to the results of 
calculations using the sudden approximation model and a planar cutoff. 



34 



1000 




2000 



1500 



1000 



500 




4800 5200 5400 

P^ [MeV/c] 



5600 



FIG. 6. A comparison of the core-fragment longitudinal momentum distributions calculated 
using the sudden model approximation with a planar cutoff (solid line) and the Glauber model [13] 
(dash-dotted line). The stripping and diffraction plus Coulomb components (thin solid and dashed 
lines respectively) in the sudden approximation are also shown. 
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FIG. 7. Comparison of the transverse (p^;) momentum distributions calculated within the sud- 
den approximation (solid line) with the Glauber model of ref. [13] (dashed line). The calculated 
distributions have been convoluted with the experimental effects. 
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FIG. 8. Momentum distribution in the projectile rest frame (panels a and c) and in the lab- 
oratory frame (panels b and d) for the s-wave ground state in ^^C (5'„ = 1.2 MeV). Top panels: 
momentum distributions in the transverse x direction. Bottom panels: momentum distributions in 
the y direction. Calculations including Coulomb dissociation are shown by the solid lines. Results 
without Coulomb interaction are represented by dashed lines. 
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FIG. 9. Intrinsic momentum distributions (dashed lines) compared with the distributions from 
the reaction calculation for an s-wave valence neutron in ^^C (top panels) and a d-wave in ^^N 
(bottom panels). Projections in the transverse x direction are shown in panels (a) and (c), whilst 
those in longitudinal direction {z) are shown in (b) and (d). 



38 



800 



600 



400 
C 200 

3 
O 

^ 

>«2500 

t 

O 2000 



1500 



1000 



500 







10000 



8000 



6000 



4000 



2000 




600 



400 



200 



-200 



200 -200 

P^ [MeV/c] 



200 



FIG. 10. Comparison of the transverse (p^) momentum distributions calculated within the 
sudden approximation with the data of ref. [13] for ^^B, 15,16q g^j^^ i9]\j_ rpj^g decomposition of the 
total distribution (dashed line) into stripping (dashed-dotted line) and diffraction plus Coulomb 
(dotted line) components is shown. The calculated distributions, after Monte Carlo filtering to 
account for the experimental effects, are also displayed (solid line). 
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FIG. 11. ^^C core-fragment transverse momentum distributions from reactions of ^^C on a car- 
bon target (data ref. [13]) calculated in the sudden approximation for three spin-parity assignments 
J^ for the ground state of ^^C. 
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FIG. 12. The same as in Fig. 10 but for the perpendicular (p_l) momentum distributions calcu- 
lated within the sudden approximation. 
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FIG. 13. Comparison of the measured and predicted longitudinal momentum distributions 
for single-neutron removal from ^^C on a carbon target. The results obtained using the sudden 
approximation model are shown by the solid line, whilst those derived using the TCM (see text) 
are displayed by the dashed line. 
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FIG. 14. Evolution of the widths (FWHM) of the longitudinal momentum distributions with 
the spectrometer angular acceptance. The experimentally observed trend (filled circles, solid line) 
is compared to the calculations in the sudden approximation model (open circles, dashed lines). 
Note that the scale for the widths has been expanded. 
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FIG. 15. Evolution of the core-fragment longitudinal momentum distribution for reactions of 
^^C on carbon as a function of angular acceptances of the spectrometer. The model calculations 
(solid lines) have been normalized to the maximum number of counts and include, in addition to 
the finite acceptances, all experimental broadening effects. 
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FIG. 16. The same as in Fig. 15 but for ^N 
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